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The notion of geometric phase has been recently introduced to analyze the quantum phase transi- 
tions of many-body systems from the geometrical perspective. In this work, we study the geometric 
phase of the ground state for an inhomogeneous period-two anisotropic XY model in a transverse 
field. This model encompasses a group of familiar spin models as its special cases and shows a 
richer critical behavior. The exact solution is obtained by mapping on a fermionic system through 
the Jordan- Wigner transformation and constructing the relevant canonical transformation to realize 
the diagonalization of the Hamiltonian coupled in the fc-space. The results show that there may 
exist more than one quantum phase transition point at some parameter regions and these transition 
points correspond to the divergence or extremum properties of the Berry curvature. 
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I. INTRODUCTION 

Since the existence of the adiabatic geometricphase 
that was first revealed in pioneer work of Berry [1| , the 
concepts of geometric phases have been extensively gen- 
eralized along many directions 2-4], and now the ap- 
plications of geometric phases can be found in various 
physical fields 5 9] . Recently, the close relation between 
geometric phases and quant um p hase transitions (QPTs) 
has been revealed gradually |10l - fl2l ] and increasing inter- 
est has been drawn to the role of geometric phases in de- 
tecting QPTs for various many-body systems [l3| . Essen- 
tially, quantum phase transitions happened at zero tem- 
perature are characterized by the dramatic changes in the 
ground-state properties of a many-body system. Unlike 
classical phase transitions driven by thermal fluctuations, 
QPTs are driven by pure quantum fluctuations. Tradi- 
tionally, QPTs are analyzed by resorting to notions such 
as the order parameter and symmetry breaking within 
the Landau- Ginzburg paradigm, however, this scheme 
cannot give a complete description of properties of the 
ground state in many-body systems. In the past few 
years, a lot of efforts have been devoted to the study of 
QPTs from other perspectives, such as quantum order or 
topological order [l4j], entanglement measures [H, Il6j . 
and quantum fidelity based on the concept of quantum 
information flrl 21) . Generally, in the vicinity of QPTs, 
the changes in the ground state driven by external param- 
eters of the Hamiltonian will lead to an energy-level cross- 
ing or avoided energy-level crossing between the ground 
state and the excited state [2l[ , and the features of level 
structures can be captured by the geometric phase of the 
ground state because the features of energy-levels cross- 
ing or avoided crossing correspond to the divergence or 
extremum property of the Berry curvature. From the 
geometrical perspective, the geometric phase is a reflec- 
tion of the global curvature in the parameter space of the 
Hamiltonian. 

In the present work, we shall use the geometric phase 
of ground state to detect the QPTs for a period-two inho- 



mogeneous anisotropic XY spin-^ chain in a transverse 
field, in which the nearest-neighbor interactions and the 
degree of anisotropy will take alternating parameters be- 
tween the neighbor sites [22T4271 ] . To the best of our 
knowledge, the previous studies of geometric phase as 
a witness of QPTs mainly concerned homogeneous spin 
chains for simplicity |lCH12j . On the other hand, inhomo- 
geneous systems will exhibit rich phase diagrams and it 
would be interesting to investigate whether the geometric 
phase is able to characterize the quantum phase transi- 
tion in these more complicate systems. So far, many 
methods have been introduced to investigate the inho- 
mogeneous spin chains in different limited conditions (22r - 
[28l |: however an explicit expression of the ground state, 
which is necessary for the derivation of the geometric 
phase, is still lacking. In our scheme, by mapping the spin 
Hamiltonian on a fermionic system through the Jordan- 
Wigner transformation and the Fourier transformation, 
we derive a general canonical transformation to realize 
the diagonalization of the fermionic system Hamiltonian 
coupled in the k- space and construct the exact expres- 
sion of the ground state. Our results show that there 
exist more than one critical points at some parameter 
region and the critical points correspond to the diver- 
gence or extremum property of the Berry curvature of 
the Hamiltonian parameter space. 



II. MODEL 



The system under consideration is an inhomogeneous 
periodic anisotropic XY spin-i chain [22T - [27l | . which con- 
sists of N cells with two sites in each cell, and in an 
external magnetic field. Its Hamiltonian is given by 

N 

JP= - Y,^ 1+ ^ S taS?, b + Ja(l-Ja)Sl a Sl b 
1=1 

+ J b (l + j b )Sf !b Sf +lta + J b (l - i b )Sf b Sf +l a 

+ h(sr a + sr b )}, (i) 
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However, the differences between the two boundary con- 
ditions are negligible in the thermodynamic limit where 
the second-order QPTs occur [TJ [29[ • Without loss of 
generality, we assume the periodic boundary condition 
on the fermionic system, which means that Np is always 
even and C7jv+i, m = Ci.m- This periodic boundary con- 
dition enables us to introduce a Fourier transformation, 



e ikR '« a k , 



K 



(4) 



where Sf m {a = x, y, z; m — a, b) are the local spin op- 
erators, J m is the exchange coupling, 7 m is the anisotropy 
in the in-plane interaction on the m site in the Zth cell 
and h is the external field. We assume periodic bound- 
ary conditions and choose TV to be odd for convenience. 
This model encompasses a group of other well-known 
spin models as its special cases 29J, such as quantum 
Ising model in a transverse field for "f a = 7& = 1 and 
J a = Jb, the transverse field XX model for j a = j b = 
and J a = J bl and the uniform transverse field anisotropic 
XY model for j a — j b and J a = Jj,. 

In order to obtain the geometric phase in this sys- 
tem, we consider rotating this model by applying a rota- 
tion of angle ip about the z axis to each spin [ldj . i.e., 

@z W) = rij^exP^C^a + S l,b% and tnen we nave to the Hamiltonian H{<p), in which k = (2tt /2aN)n and 
Jtffa) = 3>t {ip),??® z {tp), in which <2> z (tp) is the rele- 
vant rotation operator and we have set h = 1 for simplic- 
ity. It can be verified that Jf?(0) = J^(tt), and J^iip) 
is 7r periodic in tp because the quadratic form about the 
x and y axes appears symmetric in Eq. ([TJ. Consid- 
ering the unitarity of the rotation operator U (ip) , the 
critical behavior and energy spectrum of the family of 
Hamiltonians parametrized by tp are obviously ip inde- 
pendent. The spin Hamiltonian can be mapped exactly 
on a spinless fermion model through the Jordan- Wigner 
transformation 



N-l 
2 



N-l 
2 



1, 



N- 



i. Here R la (R w = i?; a + a) 
is defined as the coordinate of site a (b) on the Ith cell 
in the one-dimensional lattice with the lattice parameter 
2a. Hence, the Hamiltonian W{<p) transformed into the 
momentum space is given by 



+ b{b k - 1)+ [(^e ika + ^e- ika )a{b k 



+ H.C.] - [(^e 24V+lfca - ^ e 2l *- lka )a- k b k 



S ta = ex P 



l-l 



I' — lm' — a.b 



H.C.]} 



(5) 
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°L6 



exp 




(2) 

where S lm — Sf m ± i5f m are the spin ladder operators 
and C; and Q ;m are the fermion creation and annihila- 
tion operators. The original Hamiltonian ^(ip) is trans- 
formed into 



This Hamiltonian has a quadratic form in fermion oper- 
ators and can be exactly diagonalized. We note that the 
Hamiltonian H v can be expressed as H v — J2 k (r\M k r k + 
h) with matrix rt = (at, a,- k , bt, b^ k ) and M k is a 4 x 4 
Hermitian matrix. Therefore, we can always find a uni- 
tary transformation matrix U which can be inserted in 
the Hamiltonian as H {p = T, k ( T k u l u kM k UlU k T k + h) 
and then transform the matrix M k into a diagonal matrix 
U k M k Xj\. That is to say, the term U k T k is equivalent to 
introducing the following canonical transformation and 
define a set of new operators, i.e., 



X ^ r ,Ja 



1=1 



■ Jbib 



l2ip Cl b Cl +lM + H.C.) 



+ h{Cl a Cu a + Cl b C h b - 1)] 



(3) 



In the fermion case, the periodic boundary con- 
ditions S N+hm = Sf m Xa = x,y,z; m 
a . b) on the spin degrees of freedom imply that 

CiV+l.m = ^V[iKY^V = lY^i=a C i, m > C l',™>]Cl,rn, in 



Which (Ei'=lEm'=o C 'l',m' C 'i',»'') is J USt tlie t0tal 

fermion number Np. Thus the boundary conditions on 
the fermionic system are Cjv+i, m — e mN * 'Ci jm , and the 
fermionic system will obey periodic or antiperiodic con- 
ditions depending on whether Np is even or odd [30| . 



7fc 



1 

3 V 



—.(e 2 ^ cos^afc + e tSk e- l ' Tk sin^al 



2 2 _fe; 



+ e ^ S in °±b k + e-^e-^e*** cos ^6+ A 
2 2 

1 (e 2 ^ cos %a fe - e^e"^ sin %aL fc 



2 
2 



e 2lv e l5fe cos ^b k + e la " sin ^6 f _ fe ) , 



_L( e -M* e io» sin %a fc + e^ 2 ^ cos fe 
v/2 2 2 fc 



2CTi- * 7 

e fc sm — o fc - e 



cos— 6L fc ) ,(6) 
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where 



-ika\ 



S k = arg ( J a e lKa + J h e 

a k = arg ( J ala e* ka - J hlh e~ ika ) 

Ck = \J J 2 + Jl + 2J a J b cos 2ka 

& = 



note that the energy spectrums A^fc > 0, A„fc > and 
A 7 fc < 0, A M fe < 0. The ground state, denoted as \g(<p)), 
corresponds to the state with the lowest energy, which 
consists of state with no rj and v fermions occupied but 
with 7 and \x fermions occupied. Explicitly, the ground 
state can be constructed as follows 



Jill + Jbll ~ MaJblalb cos2fca , (7) \ g ((p)) = C~* JJ (jl k 7t^_ k ljiv-kVk^-k^k) |0) a ® |0) b , 



and 



cos e/fe = 



cos/3 fc = 



(/i-%) 2 + (%) 2 



h + % 



(/^ + %) 2 + (%) 2 



(8) 



Using the set of quasiparticle operators jkiVkt U>k and v k , 
we can write the Hamiltonian H v in the explicit diagonal 
form as 



9=7. 1.M.f fe 



(9) 



where A gi fc (g = j,rj,iu,,v) are the eigenvalues of the 
Hamiltonian matrix Mfc, and now, they are the relevant 
quasiparticle energy spectrums as follows: 
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-(f) 2 




2 7 




j(h + 


y) 2 - 


-(f) 2 



A„fc = 



Furthermore, it can be verified that the general canoni- 
cal transformation Eq. © can be reduced to the famil- 
iar Bogoliubov transformation in the case of the uniform 
anisotropic XY model. 



III. GEOMETRIC PHASE AND QUANTUM 
PHASE TRANSITION 



k>0 

(11) 

where C 1//2 is the normalized factor and |0) a and |0)& are 
the vacuum states of the sublattices a and b, respectively. 
It is easy to check that rjk\g(<p)) — 0, Vf.\g(<p)) = and 
7jtlff(v)) — 0, [i\.\g(ip)) = for all k. The corresponding 
ground-state energy E g is 



E, 



= E( A 



A 



jj.k 



h) 



(12) 



The geometric phase of the ground state, denoted 3B g , is 
given by 



;— 

dip 



g (if) ) dip 



(13) 



Substituting Eq. (fTTjl into Eq. ([13]) , we have 

1 -r-r / 

J ° k>0 



8 



i ^]l{l-klll 1 -kt J 'l r )-kVkV-kVk) |0) a ® |0) 6 dp.(14) 

^ fc>0 

The factor of | before the normalized factor C _1 is due 
to the repeated calculations about the k and —k oper- 
ators. Straightforward calculation is tedious. Neverthe- 
less the result can be derived concisely from the follow- 
ing consideration. We note that for each term of 7fc|p7| 
and 7-fc^7lfc in the integrand yield the same results of 
— 2icos 2 In the same way, the terms of Mfc^A 1 ! and 
/i-fe^/il fe yield the results of — 2icos 2 the terms of 
Vk^Vk and r] ] _ k -^r]_ k yield the results of -2icos 2 ^, 
and the terms of v\-j^v k and ^ k -§^^-k yield the results 
of — 2i cos 2 4r- Finally, the overall result is 



Now, let us focus on the geometric phase of the ground 
state. We have introduced the family of Hamiltonians pa- 
rameterized by ip, and this family of Hamiltonians H(ip) 
can be described as a result of adiabatic rotation of the 
physical system. The geometric phase of the ground state 
will be accumulated when the system finish a cyclic evo- 
lution, corresponding to varying the angle ip from to it 
[Jf?(ip) is 7r periodic in ip]. 

The Hamiltonian H v in Eq. ([§]) has been diagonalized 
in the set of quasiparticle number operators, which allows 
us to determine all the eigenvalues and eigenvectors. We 



i 
2 

2?r 



fc>0 



2 °k 



2 Pk x , 

i cos — ) dp 



(A^-l) + ^(cos6» fc +cos/3fc) 

fe>0 

27r ^2 ( cos 8k + cos (i k ) . 

k>0 



(15) 



To study the quantum criticality, we are interested in 
the properties under the thermodynamic limit when the 
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FIG. 1: (color online) (a) The derivatives of the geometric 
phase d/ig/dh for an inhomogeneous periodic Ising model in a 
transverse field h, as a function of the Hamiltonian parameters 
a = Jt/Ja and h, in which (J a = 1,7a = Ti = 1). The curves 
correspond to different lattice sizes TV = 51; (b) N = 101; (c) 
N = 501; (d) N -> oo. 

size of the spin lattice N — > oo. In this case, we in- 
troduce the notation of the geometric phase density as 
/3 g — liniAr^oo SSg/N , thus, we have 

27]- 

P g = lim — y (cos0 k + cos/3 fc ) 

TV— >oo iv ^ — ' 

fc>0 

f-TV 

= / (cos^fc + cos/3fc) dfc . (16) 
Jo 

In this case, the summation ^ J2k>o nas been replaced 
by the integral Jq dk with dfc = limAr^oo To bet- 
ter understand the QPTs of this inhomogeneous periodic 
model and how the geometric phase of the ground state 
is used as a witness to detect them, we present numerical 
results for the derivative of its geometric phase dh(3 g as 
a function of different parameters ( J a , Jb,Ja,Jb) in the 
Hamiltonian. 

In Fig.l, we plot it as a function of a = J / J a and h 
with fixed parameters J a = 1 and 7a = 76 = 1, which 
describes an inhomogeneous periodic Ising model in a 
transverse field h. As shown in Fig.l, the peak of curves 
for dhPg{a, h) becomes sharp with the increasing of the 
lattice size N. A notable feature is that the divergence 
of the curve in the thermodynamic limit only exists in 
the parameter region of Jb/J a — 1 and 7 a =75 = 1, 
which correspond to the case of the uniform quantum 
Ising model, while in the other parameters regions, the 
curves only show extremum points. 

In Fig. 2, we illustrate the derivative of the geometric 
phase of the ground state in various cases of inhomoge- 
neous periodic systems. An interesting thing is that there 
may exist two quantum phase transition points in some 
parameter regions for the inhomogeneous period-two sys- 



tems. The number of transition points and the corre- 
sponding divergence or extremum properties of curves 
are dependent on the parameters of the Hamiltonian, 
which is quite different from those of the quantum Ising 
model and anisotropic XY model in a transverse field 
[Tol rTlj |. As shown in Figs. 2(a) and 2(c), the deriva- 
tives of the geometric phase only display the extremum 
instead of the divergence properties even under the ther- 
modynamic limit condition. On the other hand, in Figs. 
2(b) and 2(d), the extremum and divergence properties 
can coexist in some parameter regions. In order to fur- 
ther understand the divergence or extremum property 
of dh(3 g (a 7 h), we choose a section of Fig. 2(c) plotted in 
Fig. 3, in which the Hamiltonian parameter takes J a = 1, 
Jb = 2, -f a = 0.05 and j — 0.08. In this case, the tran- 
sition point of QPT is characterized by the extremum 
point. 

As shown in Fig. 3, there is no real divergence even 
in the thermodynamic limit but it tends to two ex- 
tremum points with the increasing of the lattice size 
N. The transition points in the thermodynamic limit 
can also be obtained by the finite-size analysis of po- 
sitions of extremum points for different size systems. 
Our results show that the position of the first extremum 
point approaches the first QPT point h c i in a way of 
h m = h c i(l — constN~ lm4 ) with the transition point 
h cl = 0.559908 and the second one approaches as h m = 
h c2 (l — const TV -1 ' 017 ) with the transition point h C 2 — 
1.47561. 




FIG. 2: (color online) (a) The derivatives of the geometric 
phase dp g /dh as a function of 7 = — 76 and h with the 
fixed parameters J a = 1, Jb — 0.5 and lattice sizes TV — > 00; 
(b) d/3 g /dh as a function of Jb and h with the fixed parameters 
J a = 1, 7a = 0.2, 7(, = 0.4; (c) d/3 s /dh as a function of 76 and 
h with the fixed parameters J a = 1, Jb — 2, 7 a = 0.05; (d) 
d/3g/dh as a function of Jb and h with the fixed parameters 

Ja = 1, 7a = 0.2, 7b = 0.1. 
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ln(N) ln(N) 



which is generated by a unitary transformation 
H(R,X) = <^(R)ff(A)^ t (R). Here, <^(R) is unitary 
and satisfies (R) , H (X)] 7^ to ensure the nontrivial 
transformation. Obviously, such a transformation keeps 
the energy-level structures invariant and the critical be- 
havior of the system is thus R independent. The eigen- 
values are only characterized by the parameter A. On the 
other hand, we note that the geometric phase adiabati- 
cally undergoing a closed path Cra in the R space can 
be expressed as 

Pg(o an x) = - [[ n g (R,x)-ds, (20) 

where il g (R, A) is the Berry curvature given by 



FIG. 3: (color online) (a) The derivatives df3 g /dh for the 
inhomogeneous periodic XY model(J a = l,Jt = 2, y a = 
0.05, 7b = 0.08) as a function of the transverse field 
h. The curves correspond to different lattice sizes N = 
51, 101, 501, oo. (b) and (c) The positions of the first ex- 
tremum point changes and tends as jV -1 004 towards the first 
QPT point h — 0.559908; The positions of the second ex- 
tremum changes and tends as TV" 1 ' 017 towards the second 
QPT point h = 1.47561. 



IV. DISCUSSIONS AND CONCLUSIONS 

As shown above, the geometry phase can be used as a 
detector for the more complicated QPTs in the inhomo- 
geneous system. This is because there exists an intrinsic 
connection between the geometric phase and the energy- 
level structure. Furthermore, similar connection is also 
reflected in the fidelity. The relation between the fidelity 
and Berry phase has been unveiled in terms of geomet- 
ric tensors 18]. The intrinsic relationship between the 
fidelity and the characterization of QPTs has also been 
studied in Ref. For a general Hamiltonian of the 

quantum many-body system undergoing QPTs given by 



Q g (R, X)=Irn 



(.9R,A|V R i/|n R , A ) x (nR,A| Vr-H'Is'r.a) 



H(X) = H a + XH X , 



(17) 



where Hi is supposed to be the driving term with A as the 
control parameter, the second derivative of the ground- 
state energy can be expressed as [l9| 



8 2 



^ 2|<n(A)|g 1 |g(A))r 



E g {X)-E n {X) 



(18) 



Here E g is the ground-state energy, n(X) are the eigen- 
states ofH(X), and 5(A) is the ground state. On the other 
hand, the geometric phase of the system can be obtained 
by introducing another parameter R to the Hamiltonian 
Eq. (H2l), i.e., 



(E n (X)-E g (X)) 2 



H(R, X) = H (R) + Aifi(R) 



(19) 



n^g V^nVV-^gy 

(21) 

From the expressions of Eqs. (fT8|) and (f2~Tj) . it is not 
hard to find that for both of them the singularities may 
come from the vanishing energy gap in the thermody- 
namic limit. In the inhomogeneous XY model, we find 
that a gapless excitation occurs only when A n k — > or 
h-vk 0, which demands £/. — > 0. Hence, this condition 
can be achieved only in the thermodynamic limit N — > 00 
and for the appropriate parameters of the Hamiltonian, 
i.e., Jala — Jblb- Apart from these special cases, there 
exist no solution for A^fc = or A u k — and a non-zero 
energy gap opened. Consequently, the Berry curvature in 
the thermodynamic limit only develops extremum points 
instead of divergence. 

In summary, we present an exact diagonalization ap- 
proach for an inhomogeneous periodic anisotropic XY 
model in a transverse field. By introducing a general 
canonical transformation, we construct an explicit ex- 
pression for the ground state, and based on this, we study 
the geometric phase of the ground state and QPTs for 
this model. Different from the Ising chain and anisotropic 
XY chain in a transverse field, the inhomogeneous peri- 
odic spin chains exhibit a richer behavior of QPTs. Our 
results show that there may exist more than one phase 
transition point at some parameter regions. In the lan- 
guage of geometric phase, detecting the QPTs of a many- 
body system driven by the external parameter A is equiv- 
alent to finding a path Cr.a in the parameter space of 
the Hamiltonian, in which the Berry curvature comes to 
the divergence or extremum points[31]. 
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